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I comment on the relation between two sampling methods for absorbing state models. It is shown 
that a certain ensemble without external field conditional to activity coincides with the unconditional 
ensemble for sufficiently small external field. The actual physical processes involved are identical 
and the derivation of the identity of the scaling behaviour relies on a single (established) scaling law. 
While the conditional ensemble by construction does not contain information about the system with 
large external field, it contains all information about the limit of vanishing external field and about 
the vicinity of the critical point: Finite size scaling as well as critical scaling in the temperature-like 
variable or in (small) external field. 



Absorbing state (AS) models suffer from the problem that the activity necessarily ceases in any finite system in 
the stationarity state [1]. In order to see a phase transition, one cannot naively probe the stationary state. Two 
approaches are widely accepted to overcome this problem: Either an external field is applied that creates spontaneous 
activity [2] or an ensemble conditional to activity is considered [3] (or technically more sophisticated [4]). Both 
methods effectively do away with the absorbing state altogether. The additional external field might be physically 
more appealing because it resolves the problem by introducing an external driving that appears naturally in the 
corresponding field theories. On the other hand, the external field represents a parameter that needs to be tuned 
in addition to the temperature-like variable that drives the transition, because certain universal behaviour can be 
obtained only for sufficiently small field [5]. The conditional ensemble has one parameter less, but requires the selection 
of "good samples" which might appear unphysical. 

The different techniques have created a rift that goes through the AS literature: Some authors use strong words 
to reject one and support the other method. In the following it is shown, however, that the external field ensemble 
coincides with a conditional ensemble. The derivation applies to all models where the external field triggers activity 
in the same way as an initial seed in the conditional ensemble. 

None of the two ensembles is superior to the other, provided they are implemented appropriately. In fact, one can 
use the same implementation and simply derive one from the other. Although the relation between the two ensembles 
appears to be rather trivial, it has significant implications in particular for the external field ensemble. Yet, so far, it 
does not seem to have entered into the AS literature. 



In the following, the relation between the moments of the activity calculated in the different ensembles is derived. 
The precise definition of the activity depends on the particular model, yet, all that matters is that the system has an 
order parameter, the activity, the instantaneous value of which vanishes for good if the system is "left unattended" , 
that is, the system hits an absorbing state after some time. An activation mechanism can trigger a spell of activity (an 
avalanche) and these "seeds" are either implemented as an external field h or as the initialisation in the conditional 
ensemble. 

Firstly, I introduce the nth moment of the activity in the external field ensemble, This is derived from a 

record of instantaneous activity pnit), see Fig. 1, in the obvious way. 



I. INTRODUCTION 



II. DERIVATION 




(1) 



where the time Tq is the time span of the entire observation. 
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FIG. 1: Cartoon of a continuous record of activity in a system with a small external field h which triggers new avalanches 
typically with some time of complete quiescence between them. The total time span covered is To = 300. The record could 
have been produced by compiling the three instances shown in Fig. 2. 




FIG. 2: Cartoon of three avalanches produced by placing "seeds" in an otherwise inactive system. The three instances could 
have been taken from the continuous record in Fig. 1. 
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Secondly, I define the nth moment of the activity in the conditional ensemble as the average taken conditional to 
the system being active at all. Fig. 2 shows three instances of conditional activity. To define it formally, pc{t;m) is 
the mth instance {i.e. the mth run) of a spell of activity in the system at time t after initialisation ai t = 0. If is 
the time the system is active, the moments are 

\P )c = M . • (2) 

Two technical remarks: Firstly, one should formally distinguish between the estimate of the average from a finite 
sample and the ensemble average. For brevity, this is not done here. Secondly, (p")c is often estimated slightly 
differently to allow for a transient [6] : One defines 

which is a temporal average from t = T to t = oo oi the average (across instances) instantaneous activity conditional 
to activity, weighted by the number of systems active at time t. A transient can then be established as the time 
To after which the average {p^^)c(T) has reached a (quasi-)stationary state (see discussion below). Transients are 
commonly used in equilibrium statistical mechanics, to suppress the effect of the initialisation. In the following, no 
transient is discriminated, i.e. {p"')c = {p^)c{T = 0) is used. 

Both, {p^^)h and {p^)c depend on the system size L and {p'^)h is, in addition, a function of the field h. However, 
the records pc (t) can be derived from pH (t) and vice versa, provided the field h is small enough. This is illustrated in 
Fig. 1 and Fig. 2: The three records in the latter Figure are derived from the former, by restarting the clock t after 
the rnth spell of activity is triggered in PH(t). On the other hand, given an ensemble of avalanches pc{t; rn), a record 
Pnit) can be constructed by introducing gaps of vanishing activity between spells of activity taken from pc{t;m). 
Without knowing whether the pc{t;m) are independent for different m, they have to be compiled sequentially. This 
property, however, will not be used in the following. 

To construct the basic observablcs pnit) from pc{t;m) correctly, one would also need to know the distribution 
of waiting times. This can, in principle, depend on the configuration of the system and might therefore vary if the 
system has more than one absorbing state. However, in the following only the average waiting time will enter. 

If the two observables are derived from each other in this way, one has by construction 

/•To M 

{pIh - T-^ / dtpUt) = T-^ / '^^Pc^^-^ = T " ^P^)^ ■ (4) 

For a large number of avalanches M the pre-factor Tq^ Y^m=\ converges to {Wm)^ {tm)c where {wm)H is the 
average time that passes between two avalanches being triggered in the external field ensemble and {tm)c is the 
average duration of an avalanche. 

In order to derive {wm) h given h one needs to know how the external field is implemented. Provided it operates as 
a Poisson process throughout the system, one has (wm)i? = where V is the volume of the system and h is the 

field (density) or fiux. It is also clear that the mechanism to trigger avalanches used in the conditional ensemble must 
correspond to the external field operating in the system with field. This scenario is the case for the vast majority 
of models considered in absorbing state phase transitions. This is particularly obvious if only one absorbing state 
is present, so that single "seeds", which arrive independently, uniformly and with constant rate everywhere in the 
system, trigger avalanches. 

Accepting the Poissonian nature of the external field, one arrives at the central result 

{p^HiSp, h; L) = hV{L) {tm)c{Sp; L) {p^ciSp; L) , (5) 

valid asymptotically for sufficiently small field. In Eq. (5) all dependences of the individual observables on the 
parameters have been stated explicitly. The linearity in small field of all momenta {p'^)h{5p, h;L)(xh is consistent 
with the literature, e.g. [5] but see [13]. The parameter 6p is the temperature-like variable that drives a transition at 
dp = 0, subject of the following considerations. 

III. DISCUSSION 

The above construction requires the field h to be so small that no new seeds arrive while an avalanche is running, 
or, more specifically, that the statistical weight of these events is negligible. Otherwise (p")c derived from (p")i/ 
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would depend on the external field, which has been excluded in Eq. (5) and is crucial for the following derivations. 
In other words, the derivation relies on a separation of time scales; in this sense (/o")c if' a-n observable as would be 
used in Self-Organised Criticality [7], while {p"')h represents the "instantaneous activity" [8, 9] 

Therefore, Eq. (5) can only be expected to hold if {Vh)~^ is large compared to {tm)c- This is the case for sufficiently 
small field at given system size, i.e. in particular for finite size scaling where the system size is fixed before the limit 
/i — > is taken. The condition is also met in the inactive phase, 6p < 0, where the average duration of the avalanches 
converges with increasing system size, so that, again, the average duration can be chosen to be small compared to the 
time scale set by the external field. In this case, the thermodynamic limit can be taken before /i — > 0. However, in 
the active phase, the average duration diverges extremely fast with system size, so that the external field can never 
be sufficiently small when the thermodynamic limit is taken first. 

Because the external field h appears only as a pre- factor on the RHS of Eq. (5), it is very simple to derive the 
scaling of {p")h{Sp, h; L) from the scaling of {tm)c{5'P\ L) and {p"')c{Sp; L) assuming, naturally, V{L) = ayL'^, where 
d is the spatial dimension of the system. The converse, i.e. the derivation of the scaling of {p")c{Sp;L) from the 
scaling of {p^)h{Sp, h; L), is not as simple and therefore shown explicitly in the following. 

Standard scaling assumptions suggest [5] 

h; L) = an\-^"R^ {a^SpX, a^hX^ aLLX-"^) (6) 

with non- universal metric factors [10] an,ap etc., arbitrary (scaling) parameter A, scaling function i?^ and the usual 
critical exponents a (field exponent), i/^ (spatial correlation length exponent) and the latter being /3i = /? (order 
parameter exponent) for the first moment and /?„ = n/3 if gap-scaling applies [11]. Comparing Eq. (6) to Eq. (5) 
implies that is asymptotically linear in the (small) external field 

(apSp, tthh, UlL) = UhhRniapSp, a^L) , (7) 

where Rn must obey Eq. (6) as well, i.e. Rn{apSp,aLL) = X'^~^^ Rn{ap5pX,aLLX~''^). 

Eq. (7) implies that the susceptibility xh{Sp, h: L) = , or the derivative of any moment (p'^)h with respect to 
h, is independent from the external field under the conditions stated earlier. 

To derive the scaling of (p")c(<5p; L), one needs the scaling of {tm)c{Sp,L) which is well understood [5], 

{tm)c{Sp,L) = atT{apSp,aLL) = atX^^^^^-^^TiapSp^aLLX-"^) (8) 

where i^n is the temporal correlation length exponent, 6 is the survival exponent and {1 — 6) = a — di>±, which will 
be crucial for the next step. Combining Eqns. (5), (6), (7) and (8) produces the scaling the scaling of the conditional 
moments: 

{pnc{6p;L) = ^h^x-''-R^iapSpX,aLLX-'^-) , (9) 
ay at 

where R^{x,y) = i?„(x, 2/)/(y'^T(x, y)). 

Comparing Eq. (9) to Eq. (6) shows that both ensembles display the same scaling behaviour. The lack of this 
proof might explain the critical remarks about the conditional ensemble in [2]. The finite size scaling for the two 
observables, however, differs: Setting 6p = and choosing A so that aLLX~"^ = 1 Eq. (9) gives 

(p")c(0; L) = ''"^'^\ -^"/-"<(0, 1) . (10) 

ayClt 

Because (p")c(0;L) is non-zero by definition but bound from above, i?^(0, 1) must be finite, so that (p")c(0;i) oc 
£j-0n/i'± a,g one would expect from equilibrium critical phenomena. Similarly, one finds 

(p")H(0,ft;L) = anal^-/''^L-f'-/-^R^{0,ahhiaLLr/-\l) 

= anal^-'''^L-f^-'''^ahKaLLy'''^Rn{Q, 1) (11) 

from Eq. (6) and Eq. (7). Again, the activity {p") h{0, h; L) is bound from above and does not vanish at finite 
external field, so that .R„(0, 1) is finite, and therefore (p")jj(0, /i; L) oc . This rather unusual finite size 

scaling behaviour of the activity is documented in the literature [5] (in particular Figures 33 and 47). The effective 
finite size scaling exponent for n = 1 is (cr--/?) = which is the finite size scaling exponent of the susceptibility. 

Indeed, for small external fields at i5p = and in finite systems, one expects {p)h = hxH- 
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Eq. (5) implies that, under the general conditions stated earlier, all moments vanish linearly in the external field. 
The Binder cumulant 

Q{Sp,h,L) = l-^^ (12) 

therefore diverges like h~^. In [2] it has been suggested that a divergent Binder cumulant is a "characteristic feature 
of all absorbing phase transitions" , but in the light of the above derivation, it appears rather like a generic feature of 
the particular ensemble. Moreover, replacing by the conditional averages {p^)c in Eq. (12) renders the Binder 

cumulant a universal moment ratio again and reinstates it as a "very useful method" to identify the transition. The 
same applies to lower moment ratios, which can be determined with higher accuracy, such as (p)^;/ (p^)c- 
It is now also clear why the slightly more complicated moment ratio proposed in [12] 



{p')h{p')h-{p)h{pYh 
{p)h{p')h-{p)h{p'')% 



TT — \^ /" ^' \^ /H /-,o\ 

^ — ~~rr\ ~\ ~ / ^2\2 y'^'^) 



converges to a finite value: Both, denominator as well as numerator are, to leading order, quadratic in h. In fact 

h^o h^o {p)h{p^)h {p)c{p^)c 

where the first equality holds because {p)h {p^)h G C'(/i^), and the second because the remaining terms carry the 
same prefactor, namely {hV{L) {tm}c{Sp- L))"^ , see Eq. (5). The findings above are also confirmed by the exact results 
in [12] (the RHS of Eq. (21) in [12] should read 2/{-kx) for x 0). 

The linearity of {p)h in h, Eq. (7), contradicts the mean field result {p)h{Sp = 0,h) = \fh [5], which is, however, 
consistent with mean field theory not allowing a proper distinction between conditional ensemble and external field 
ensemble. 

The limitations of Eq. (5) become clearer when comparing to the standard scaling assumption limi^oo {p) h{5p = 
0, h; L) DC h^/'^ [1]. Eq. (5) seems to suggest that (3/a = 1 which clearly is not the case. However, Eq. (5) appHes only 
where the waiting time [hV)^^ oc L~'^ is large compared to {tm)c oc i'^/'^-i--'^, which cannot be the case if a/v±_ > 
and the thermodynamic limit taken at finite field h. It is worth noting, however, that Eq. (7) relies solely on the 
Poissonian nature of the external field. 

Finally, I want to discuss briefly two features of conditional ensembles as they often appear in the literature, namely 
the initial condition and the discounting of a transient. The latter was introduced as T in Eq. (3), by defining the 
conditional nth moment as an average taken from t = T to t ^ oo. For any non- vanishing transient T > the 
strict identity of conditional and external field ensemble breaks down. To restore the above results, more elaborate 
arguments are needed, based on the length of transient, its contribution to the average and its scaling. 

Similarly, there is no strict identity of the two ensembles, if the conditional ensemble is based on initial conditions 
which do not correspond to configurations reached when the external field operates on inactive configurations in the 
external field ensemble. On the other hand, one might hope that the introduction of a transient erases the dependence 
on the initial condition and that asymptotically, the effect of the transient and the initial condition becomes negligible. 



IV. SUMMARY 



To summarise, I address the issue of the relation between the two ensembles usually considered in absorbing state 
phase transitions. It turns out that the ensemble of "activity conditional to activity" can be derived from the ensemble 
obtained by applying an external field and vice versa, provided that the avalanche duration is small compared to the 
waiting time between two avalanches. Technically, both methods can be implemented in (almost) the same way. 
Because moment ratios, such as the Binder cumulant, remain finite in the conditional ensemble and because moments 
of the activity in the presence of an external field all vanish linearly in the external field, the conditional ensemble 
might be more comfortable to study. 
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